Bichromatically driven double well: parametric perspective of the strong-field control 

landscape reveals the influence of chaotic states 
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The aim of this work is to understand the influence of chaotic states in control problems involving 
strong fields. Towards this end, we numerically construct and study the strong field control landscape 
of a bichromatically driven double well. A novel measure based on correlating the overlap intensities 
between Floquet states and an initial phase space coherent state with the parametric motion of the 
quasienergies is used to construct and interpret the landscape features. "Walls" of no control, 
robust under variations of the relative phase between the fields, are seen on the control landscape 
and associated with multilevel interactions involving chaotic Floquet states. 



I. INTRODUCTION 

Controlling quantum phenomena in atomic and molec- 
ular systems is an active area of research today. Persis- 
tent theoretical and experimental efforts have led to the 
formulation of several ingenious control schemes^' 3 -'^ for 
both atomic and molecular systems. Recently researchers 
have been interested in determining the "control land- 
scape" for specific systems^^ ^ 10 ] 11 ! 12 ] 13 . The moti- 
vation for undertaking such extensive studies is atleast 
twofold. Firstly, several optimal control efforts^&i^ seem 
to succeed despite employing local search algorithms in a 
prohibitively high dimensional control parameter space. 
Rabitz and coworkers^ reason that the topology of the 
control landscape can shed light on the robustness of the 
control strategies. Secondly, many control schemes 3 - in- 
volve few, carefully chosen, levels which provide valu- 
able insights into the control process. However, such 
few-level schemes might be compromised due to several 
reason o 15 ' 16 i 17 i 18 i 19 and one hopes that identifying the 
features of the control landscape associated with the de- 
viations might lead to better control techniques. 

In order to set the stage for our work we briefly discuss 
the concept of a control landscape. Consider a Hamilto- 
nian of the form H = Hq — fj£(t) with /i being the dipole 
operator and £(t) being the control held. There are two 
perspectives on the control problem which are of inter- 
est to the present work. The first one involves using a 
given field £(t\ A) characterized by one or more control 
variables A = (Ai , A2, . . . , A/). The objective in this case 
is to study the variation of, say, the transition proba- 
bility P if = \(f\U(T,0)\i)\ 2 between an initial state \i) 
and a final state |/) at some target time T. The con- 
trol landscape is constructed by viewing the transition 
probability as a functional of the field parameters i.e., 
Pif — Pif[£(t; X)]. Regions on the landscape with high 
values for Pj f yield the best possible choice of the con- 
trol variables in the presupposed £(t; A). Note that it is 
quite possible to have several distinct choices for the con- 
trol variables which yield the same value for Pi f . Several 
such examples of control landscapes can be found in the 
literature^&i&i 3 . 

The second perspective on the control problem comes 
from the optimal control theory (OCT) approach^. In 



OCT a cost (landscape) function of the form <&([/) = 
Tt[Up(0)W&\ = (9), with p(0) being the initial den- 
sity matrix and an arbitrary observable O, is extrem- 
ized with respect to the field £{t). The optimal fields 
correspond to the solutions of the variational condition 
5$>(U)/5£(t) = wherein the unitary time evolution op- 
erator U = U[£(t)] is naturally considered as a functional 
of the field. The transition probability Pif is a special 
case of $([/) with p = and 6 = |/)(/|. Again, 

the set of optimal fields that yield the same (O) consti- 
tutes a control level set whose topology is of considerable 
interest since one can imagine 'morphing' an experimen- 
tally undesirable optimal field to a more desirable field. 
In order to be able to carry out such a morphing it is 
essential to have full control of the system. Assuming 
full controllability Rabitz and coworkers have explored 
the topology of the OCT landscapes in great detail^. 
A particularly striking result is that the full quantum 
control landscape is devoid of false traps or suboptimal 
extremal. However, in practice one is invariably forced 
to impose constraints on <&(U) in various forms. Exam- 
ples include specification of target time T, minimizing the 
field fluence, and reduced or restricted parametrizations 
of the control field. 

Among the two perspectives on control outlined above 
we adopt the first one in this work i.e., the nature of 
the control landscape with a given form of the field. 
Although we do not undertake an OCT calculation in 
the present work, it is nevertheless useful to think of 
the field used herein as a result of a reduced dimen- 
sional parametrization of an optimal field coming out 
of a OCT calculation. The reason being that restricted 
parametrizations of the optimal £(t) typically projects 
the full OCT control space to a lower dimensional space 
which can exhibit suboptimal features^. 'Suboptimal' 
features on the control landscape are, therefore, expected 
in our analysis. Clearly, in the OCT context, suboptimal 
features can be avoided if the restricted parametrization 
is done carefully by preserving the overall intrinsic topol- 
ogy of the control landscape and level sets. Discovering 
such appropriate parametrizations might not be easy but 
a recent combined experimental and theoretical work by 
Roslund et ai. provides a clear example 2 °. In certain 
instances, as in the work by Bartelt et ai, a systematic 
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study of simple parametrizations of the optimal field can 
provide mechanistic insights into the control process 21 -. 

Based on the discussion above, one can interpret the 
suboptimal features as artificial due to inappropriate field 
parametrization or as genuine due to the unavoidable 
constraints imposed on the system. Consequently, the 
model study presented in this work incorporates two of 
the constraints that are expected to result in a com- 
plicated control landscape - strong fields and simple 
forms of the control field. It has been suggested that in 
strong field regimes local traps may appear more read- 
ily. However, even in the strong field regimes recent 
studies indicate that interpreting the various features of 
the constrained field control landscape is far from being 
straightforwar d 10 ! 11 ' 12 . In this context, it is well known 
that the classical phase space of strongly driven systems 
exhibits widespread chaos. Could it be that local traps on 
the strong field control landscapes arise due to the chaos? 
Unfortunately, not much is known about the impact of 
chaotic or delocalized quantum states on the topology of 
control landscapes and level sets. Such an analysis for the 
full OCT landscape is far from being simple. Therefore, 
in this work we make a beginning, within the perspective 
discussed above, to understand the nature and origins of 
the suboptimal features from a classical-quantum corre- 
spondence viewpoint. 

In most physical systems the phase space is mixed 
regular-chaotic and there is sufficient evidence in the lit- 
erature to indicate that quantum states delocalized in 
the chaotic regions can significantly influence the con- 
trol of initial states localized in the regular regions. For 
instance, Na and Reichl 1 — have shown that adiabatic pas- 
sage mechanisms are modified due to the presence of 
chaotic states leading to chaos-assisted passage. Gong 
and Brumer-1§> have recently reviewed the subject of co- 
herent control of quantum phenomenon in the presence 
of chaos. Their studies, especially using the (modified) 
kicked rotor, have highlighted several intriguing features 
like faster-than-classical quantum anomalous diffusion 
and the nontrivial influence of classical phase space struc- 
tures with area much smaller than the effective Planck 
constant. As a last example we mention the work by 
Takami and Fujisaki 2 ^, on optimal control of random ma- 
trix systems. Random matrix systems correspond to 
strongly chaotic classical phase space with no regular re- 
gions. In such cases it is possible, using OCT, to ob- 
tain an analytic expression for the optimal field which 
controls the transition between two random states. The 
aforementioned studies, and other earlier investigations, 
concerned with the control of quantum processes in clas- 
sically chaotic systems have raised many interesting is- 
sues. However, there are very few attempts to analyze 
the various features on the control landscape in terms of 
the phase space nature of the quantum states. Control- 
lability of quantum dynamics despite the strong chaos 
is surprising and calls for a careful study of the control 
landscape from a classical-quantum correspondence per- 
spective. 



We begin in section [IT] with a brief description of the 
model system, driven double well, which has served as a 
paradigm for many coherent control studies. The control- 
lability of the dynamics using a bichromatic field is stud- 
ied as a function of field parameters. The role of chaotic 
(delocalized) Floquet states in attempting to control the 
decay of an initial state localized in the regular region 
of the phase space is highlighted. Section IIIII introduces 
a measure which involves correlating the response of the 
Floquet states to the control field with the overlap in- 
tensities between the initial state and the Floquet states. 
We show that the correlator is related to generic land- 
scape functionals and can therefore yield useful insights 
into the features of the control landscape from a classical- 
quantum correspondence point of view. In section HVl we 
show that the proposed correlator is very sensitive to the 
underlying phase space structure and clearly identifies 
regions of no control involving delocalized (chaotic) Flo- 
quet states. Section fVl concludes with comments on the 
role of the chaotic states in the process of control and 
future studies. 



II. MODEL SYSTEM 

The model of choice to illustrate our approach is the 
driven double well system. The Hamiltonian is given by 

H{x,p;t) = J-p 2 + Bx 4 -Dx 2 +x£{t) (1) 

with £{t) being the field. The driven double well model is 
an ideal testbed for several schemes that have been sug- 
gested for controlling coherent superposition of quantum 
states 2 ^ 13 ' 23 ' 24 ' 25 ' 26 . Apart from being used to elucidate 
phenomena like chaos-assisted tunneling— and coherent 
destruction of tunneling 2 ^, the model above is also rel- 
evant in research areas like quantum computing 2 - and 
quantum dots 2 ^. In this context we also note an early 
study 3 - - by Schwieters and Rabitz on the above system 
with the aim of designing an optimal £{t) to reduce the 
difference between quantum observables and their clas- 
sical analogs. In the present work, as alluded to in the 
introduction, our goal is to study control using a pre- 
supposed field and the resulting influence on the control 
landscape. A variety of choices can be made for £ it) and 
in this work we choose a bichromatic form 

£(t) = Ai cos(o;<) + A2 cos(2cj£ + 4>) 

The above choice for £(t) is motivated by a large number 
of studies 3 -! that have shown the utility of such fields in 
controlling the dynamics with the possibility of using the 
relative phase 4> as an extra control knob. Moreover, due 
to the periodicity of Eq. [1] H(x,p,t) — H(x,p,t + Tf), 
with period Tf — 2tt/uj, the dynamics can be conve- 
niently analyzed 33 using Floquet theor y 34 ' 35 in terms of 
the Floquet states {|Xn)} and the associated quasiener- 
gies {E n }. 
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Throughout this work h = 1 and we fix the parameter 



value; 



,26,27 



1,D = 10, B = 0.5 and, u = 6.07. The 



double well supports about eight tunneling doublets and 
u) nearly corresponds to the energy separation between 
the ground and first excited states, neglecting the tunnel 
splittings. The primary field strength Ai € [6.0,11.0] is 
chosen such that the classical phase space, for A2 = 0, 
has two symmetry related regular islands embedded in 
a chaotic sea (cf. inset in Fig. UK). The symmetry re- 
lated regular regions and the associated symmetric pairs 
of Floquet states exist due to the discrete symmetry 
H (p, x, t) = H(—p, —x, t + n/u>). The parameter values 
have been chosen to make contact with an earlier work 
by Farrelly and Milligan 26 which explored the possibil- 
ity of controlling the tunneling using bichromatic fields. 
Different choices for the parameters, as long as the phase 
space exhibits substantial chaos, show similar behaviour 
and do not alter our conclusions. 
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FIG. 1: (Color online) (A) Decay time T versus the field 
strength Ai(A2 = 0) for the coherent state, \z) = \xo,po) = 
I — 1.5, 0), localized in the left regular island (indicated by an 
arrow in the inset showing the stroboscopic phase space at 
Ai = 6.4). (B) and (C) show the influence of the 2oj-control 
field with <f) — on the survival probabilities for Ai = 9.2 
and 10 respectively. Increasing A2 from (black) to 0.2 (red 
circles), 0.7 (blue squares), and 2.1 (green triangles) leads 
to localization in case (C) but not in case (B). The same 
information, with consistent color code, is shown in (A) as 
points unaffected by increasing A2. 



Control landscape: influence of the chaotic 
states 



The control problem that we address is as follows. The 
initial state, as in the earlier studie o 26 ' 27 , is a coherent 



state \z) = \xq,pq) localized in the phase space about 
(xo,j>o) within one of the regular islands. The coherent 
state is constructed as linear superposition 



b=0 



cb (xo,po) I b) 



(2) 



of a harmonic oscillator basis {\b}} with frequency ujb — 
6.25 (corresponding to ((Ax) 2 ) = 0.08). The coefficients 
are given by 



(3) 



c b (x QlPo ) = 4=a b e -H 2 /2 e — 0P0/2 



lu> b / I 

— [ Xo -\ po 

2 V LOb 



We employ the same basis to compute the Floquet states 
and quasienergies used in the rest of this work. The time 
evolved state \z(t)) can be expressed in terms of the Flo- 
quet states as 

\z(t)) = Y,^ Ent \Xn(t))(Xn(0)\z(0)) (4) 



Measuring time in units of Tf and owing to the period- 
icity of the Floquet states, \xn{t)) = \xn(t + Tf)), the 
above equation simplifies to 



-iE„ 



|X n (0)><Xn(0)|*(0)> 



(5) 



with t = kTf and integer k. In the absence of the 2w-field 
("control field") \z) decays over a timescale T which is 
determined from the first vanishing of the survival prob- 
ability 



S(t) = |<z(0)|z(r))| : 



-iE r 



— ^ PznPzm& 



T (z(0)\ Xn (0))( X n(0)\z(0)) 

i(E„-E m )T 



(6) 



where the overlap intensities are denoted by p zn = 
|(z(0)|xn(0))| 2 . Note that S(t) is just a special case of 
Pif with \i) = \f) = \z). In Fig. [TJA the decay time is 
shown as a function of Ai for A2 = 0. The strong fluctu- 
ations in T, despite similar phase space nature over the 
entire range of Ai, is a direct consequence of the under- 
lying chaos 2 —. Recently it has been establishe d 36 ' 37 that 
the local structure of the phase space about the location 
of \z) can play a crucial role in the decay dynamics. The 
local phase space structure around the left regular island 
is shown in Fig. [5] for various values of Ai. Comparing 
the features of the decay time plot in Fig. [TJ\ with the 
local phase space nature in Fig. [2] shows qualitative cor- 
respondence. In particular, the drastic variation of T 
around Ai = 8, seen in Fig. [TJ\, is reflected in the signif- 
icant change of the local phase space nature in Fig. [2£b) 
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convenience; the decay time is a better choice which re- 
quires considerable effort but the gross qualitative fea- 
tures of the control landscape do not change upon using 
(S). Large (small) values of (S) indicate that the decay 
dynamics is suppressed (enhanced). It is clear from Fig. [3] 
that the landscape is highly convoluted - a consequence 
of the simple bichromatic choice for £(t) and the non- 
linear nature of the classical dynamics. An earlier OCT 
study2£ on a different classically chaotic system has also 
suggested that the control landscape could be highly con- 
voluted. From a control point of view, there are regions 
on the landscape for which a monotonic increase of A2 
leads to increasing control. At the same time there are 
regions, especially around Ai » 9.2, where increasing A2 
hardly effects the decay dynamics and indicates very lit- 
tle control. As noted before, one expects such features 
even on constrained OCT landscapes. In this work we 
would like to understand the nature and origin of such 
robust regions of no control. 



FIG. 2: The local nature of the phase space in the vicinity of 
the initial state \z) (black circle) for Ai equal to (a) 7.0 (b) 8.0 
(c) 9.0 and (d) 10.0 with A2 = 0. Note the drastic change in 
(b) which is mirrored in the decay time behavior in Fig. [TJ\. 

as compared to those in Fig. [2ja),(c), and (d). There 
is also preliminary evidence (not discussed in this work) 
which indicates that the various field-matter nonlinear 
resonances control the decay of \z) and partly responsi- 
ble for some of the structure seen in Fig. [TJ^. 

Is it possible to control the decay of \z) with an appro- 
priate choice of the parameters (A 2 ,0) of the 2w-field? 
In this instance the problem of control, in the asymp- 
totic limit t — > 00, is equivalent to the suppression of 
tunneling between the two symmetry related islands in 
the phase space. At the outset one expects, for A 2 7^ 
and ^ 7r/2, the discrete symmetry to be broken leading 
to localization and hence perfect control. However, the 
results of the computations summarized in Fig. [1] show 
a more complicated behaviour. For most of this study 
we will choose two representative cases to illustrate the 
issues. In the first case, corresponding to Ai = 10.0 and 
labeled as (C) in Fig. [I] a monotonic increase of A2 with 
= results in the localization of \z). This is in ac- 
cordance with the earlier observations by Farrelly and 
Milligan 26 . On the other hand, the decay dynamics for 
the second case, corresponding to Ai — 9.2 and labeled 
as (B) in Fig. [TJ is unaffected by the control field even 
for fairly large values of A 2 . 

To obtain a global picture, in Fig. [3] the control land- 
scape for 4> = is shown in terms of the time-smoothed 
survival probability 

(S) = - f ' dr'S{r') (7) 
T Jo 

as a function of the field parameters (Ai, A 2 ). Note that 
the choice of (S) to represent the landscape is made for 




FIG. 3: (Color online) Time smoothed survival probability 
representation of the control landscape for <j) — 0. The initial 
state is the same as in Fig. [1] Notice the convoluted form of 
the landscape with the regions of low probability indicating 
little to no control. 

Insights into the lack of control can be obtained by 
studying the variation of the Floquet quasienergies with 
A 2 . In Fig. [4j\ we show the results of such a compu- 
tation corresponding to region B (cf. Fig. [T]). Over 
the entire range of A 2 six states contribute to the de- 
cay of \z) as evident from the overlap intensities shown 
in Fig. 2)3. A closer look reveals the existence of states 
delocalized in the chaotic regions of the phase space - 
two such states are shown in Fig. [4)3 in terms of their 
Husimi distributions^. In contrast, Fig. HP shows that 
the dynamics of case (C) in Fig. [T] is dictated by local- 
ized states which explains the ease of control with the 
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FIG. 4: (Color online) (A) Variation of the Floquet quasiener- 
gies with A2 and fixed Ai = 9.2. Six states that have appre- 
ciable overlap with \z) are highlighted by circles. (B) Overlap 
intensity p zn for Ai = 9.2 and A2 = 2.1 indicates multilevel 
interactions involving the states shown in (A). Husimi dis- 
tribution function of two delocalized Floquet states (red) are 
also shown. (C) Same as in (B) with Ai = 10 and A2 = 0.7 
shows the localized nature of the participating states. 



symmetry breaking 2u-field. Therefore it is reasonable 
to expect that the lack of control, signalled by plateaus 
in Fig. [TJ'V and the valleys in Fig. [3l is due to the dom- 
inant participation by chaotic states i.e., chaos- assisted 
tunneling 38 . The plateaus arise due to the fact that the 
coupling between the localized states and the delocal- 
ized states vary very little with increasing control field 
strength A2 - something that is evident from the Floquet 
level motions shown in Fig. |5J^ and established earlier by 
Tomsovic and Ullmo in their seminal work^ 8 - on chaos- 
assisted tunneling in coupled quartic oscillators. It is 
important to note that for A2 = the chaotic states, as 
opposed to the regular states, do not have a definite par- 
ity. Consequently, presence of the 2w-field does not have 
a major influence on the chaotic states^. Thus, if one or 
more chaotic states are already influencing the dynamics 
of \z) at A2 = then the bichromatic control is expected 
to be difficult. 



Similar conclusions were arrived at in a previ- 
ous study-i£ involving bichromatically driven pendulum. 
However, several questions, not addressed by the previ- 
ous studies, arise at this juncture. What is the precise 
role of the chaotic states for control? Are the plateaus 
robust for varying A2 and 0? Is there a suitable measure, 
generalizable to systems with higher degrees of freedom, 
that can identify regions on the control landscape with 
little or no control? 



III. SENSITIVE MEASURE FOR THE 
CONTROL LANDSCAPE 

In order to address the questions raised above we first 
suggest a measure, different from S(t) or (S), to map out 
the control landscape. Note that the Floquet states and 
quasienergies of Eq. [1] are parametrically dependent on 
the field strengths i.e., E n = E n (X) and \x n ) = \XnW)- 
Thus the response of a Floquet state to the 2w-field, as 
shown in Fig. [4j\, is measured by the parametric deriva- 
tive dE n /d\2 which is also called as the level "velocity". 
At the same time the decay dynamics of a coherent state 
\z) is dominated by Floquet states that have appreciable 
overlap with \z) (cf. Eq. [6]). Therefore, Floquet states 
\Xn) with substantial overlap p zn and finite dE n jd\i are 
expected to be important in controlling the dynamics of 
\z) . The qualitative argument can be made quantitative 
by introducing the overlap intensity-level velocity corre- 
lator 

C.{Xl,X) = — (pJ^®) (8) 

C Z (X) was first proposed by Tomsovio 3 ^ as a sensitive 
measure for deviations from ergodicity in strongly chaotic 
systems and has been studied in detail in subsequent 
works 40 . In Eq. [H <r z and a E are the local variances 
of p zn and dE n /d\2 respectively and the average is over 
all the Floquet states. 

The sensitivity of the correlator to phase space struc- 
tures has been amply demonstrated in the previous 
studies 40 . Since the earlier works did not explicitly deal 
with driven systems, in Fig. [5] we give one example of the 
sensitivity in the present context. In Fig. [5] we place co- 
herent states at different phase space points on the phase 
space and obtain a correlator map. Note that the coher- 
ent state parameters are identical to the ones used in 
the earlier sections. The correlator map is shown for two 
different field strengths (A) (Ai,A 2 ) = (9.2,0.7) and (B) 
(10.0, 0.7) with (j> — 0. We remind the reader that these 
parameter values correspond to the two exemplary cases 
denoted as (B) and (C) in Fig.QJA). Several observations 
can be made by inspecting the figure and we mention a 
few important ones. Firstly, the inhomogeneous nature 
of the phase spaces is reflected by the correlator. Specif- 
ically, the large values of C z in the regular islands on the 
periphery of the phase space imply participation of "fast 
moving" Floquet states in the dynamics. Examples of 
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FIG. 5: (Color online) The intensity-velocity correlator (cf. Eq. [8]) computed over the entire phase space at (A2 = 0.7, = 0) 
for (A) Ai = 9.2 (B) Ai = 10.0. Calculations are performed by placing coherent states at the different phase space points. Note 
the similarity of the phase space but the essential differences in the correlator maps. See text for details. 



such Floquet states can be seen in Fig. QJA) as states 
with large level velocities. Secondly, the correlator takes 
on different values even in the chaotic regions and there- 
fore points to the different nature of the dynamics i.e., 
different groups of Floquet states that influence the dy- 
namics. Finally, despite the similar looking phase spaces, 
it is clear from Fig. [5] that the correlator maps are dif- 
ferent. In particular, note that the correlator values in 
Fig. [SJ^ about the left and the symmetry related right 
regular islands have the same sign. However, in Fig. [5)3 
the same regions have correlators of opposite signs. This 
indicates a break in symmetry and hence provides some 
hints to the fact that one obtains control in this case as 
opposed to in case A. 

Before proceeding further we note that the correlator 
is related to the generic landscape functionals employed 
in OCT studies 6 -. In order to see this consider the generic 
quantum control landscape functional 6 

&(U) = Tr [U(t)p(0)U(t^Q] (9) 

with p(O),0, and r being the initial density matrix, an 
arbitrary observable operator, and final dynamical time 
respectively. Choosing p = \z)(z\ and = £2 U (x,t) = 
x cos(2djt + <f>) one can express the landscape functional 
as 

HU) = ^(Xn(0)|z(r))(z(r)|£ 2 J Xn (0)) (10) 

n 



Using Eq. [5] one can write 

$(U) = ^e- ! ( £ "- £ ™)^ Xn |z)(z|A m )(x m |^(r)|Xn> 

rn,n 

fdEJX)\ 

= / Pzn — t: + nondiagonal terms (11) 

' V 0X2 J 

n x ' 

The first term on the last line in Eq. Illl comes about due 
to the Hellman-Feynman theorem. The diagonal term in 
Eq. [11] is essentially the correlator in Eq. [5] and the con- 
tribution from the nondiagonal terms can be neglected 
in the long time limit. The variances in Eq. [H] ensure 
that C z is a true correlation function. However, without 
the variances one is dealing with the covariance which is 
related to the landscape functional <E> via Eq. QT] Note 
that the covariance already has the essential features that 
are needed for our study and the existing semiclassical 
insights 40 should prove valuable in the context of the 
present work. 

What should one expect from the correlator? In gen- 
eral, C Z (X) highlights the presence of groups of states ex- 
hibiting common localization feature o 39 ' 40 . In the ergodic 
limit it is known^ 9 . that C z (\) ~ ± N^ 1 / 2 for any ini- 
tial state I z) with TV being the number of states involved 
in the averaging in Eq. [H] However, despite the strong 
fields and extensive chaos, the system considered herein 
is not ergodic (cf. Fig. [5]) . The interpretation is subtle as 
compared to the earlier studie o 39 ' 40 since chaotic states 
already exist for A2 = 0. The chaotic states, due to their 
delocalized nature, have small level-velocities \dE n / 8X2] 



7 




-0.2 0.3 0.8 400 1 800 

FIG. 6: (Color online) Control landscapes constructed via the intensity-level velocity correlator C Z (A) for phase cj> equal to (A) 
(B) 7r/4 (C) 7t/3 and (D) n/2. All landscapes shown have the same axes range. Regions lacking control are associated with 
Cz(A) < 0. Note the one-to-one correspondence between the landscape in (A) and the one shown in Fig. [3] and the robustness 
of the "wall" of no control around Ai ~ 9.2 (indicated by arrows). Survival probabilities shown in the lower panels for specific 
field parameters (indicated by circles) confirm the predictions based on C 2 (A). Calculations of the landscape are performed on 
a 50 x 30 grid with a spacing of 0.1. 



which is evident from Fig. 3)A_ as well. Consequently, par- 
ticipation of several chaotic states in the decay dynamics 
of \z) results in small value for |C Z (A)|. As a simple 
model, assume iV c chaotic states with nearly the same 
level velocities v c contribute equally to the dynamics of 
1 2) i.e., similar intensities p zc — N~ . Note that the di- 
vision of the total number of states as N — N c + N r , 
with N r being the number of regular states, is strictly 
valid only in the semiclassical limit. Clearly the param- 
eters in this work are such that the system is far from 
the required limit but the arguments are presented with 
an attempt to understand some of the features of the 
correlator. Within the simple model, using the sum rule 
that the average of the level-velocities over all N states 
is zero, it is possible to show that 



c z = 



3 = 1 

/ch(l + /ch) 
«2/ch + (V 2 ), 



fc/ch 

1/2 



,2\l/2 



(12) 



(13) 



with / c h = N c /N r and the limiting form of the correlator 
above arises when v~ 2 (v 2 ) r > /d, > 1. It is crucial 
to observe that the sign of the correlator in Eq. [T3] is 
determined by the chaotic level- velocity v c . An a priori 
knowledge of the sign of v c in a specific system is difficult 
but can be determined from a single Floquet calculation. 



Eg. fT51 suggests that \C Z \ will be small for (v) r <C 



.,2\!/2 



IV. INTERPRETATION OF THE FEATURES 
ON THE CONTROL LANDSCAPE 



Armed with the qualitative insight into the behavior 
of C z {\) we now compute the control landscape using 
Eq. [5] for the entire range of Ai shown in Fig.[T]A_. The re- 
sulting control landscape is shown in Fig. [SJ'V for relative 
phase 4> = and should be compared to the landscape in 
Fig. mapped via the time-smoothed survival probabil- 
ity (S). It is important to note that the two landscapes 
show identical features. In particular, the regions lacking 
control on the landscape in Fig. [3] correspond to the re- 
gions showing C Z (A) < in Fig. OV. However, using the 
correlator based landscape it is now possible to provide 
an explanation, not readily forthcoming from the land- 
scape in Fig. [3l for the source of such lack of control. For 
the specific parameters chosen in this study the chaotic 
level-velocity distribution is peaked around a negative 
value and thus C z (\) < 0, following the results of the 
simple model in Eq. [131 implies that the regions of no 
control arise due to the participation of chaotic Floquet 
states. Although the role of chaotic states to control is 
already highlighted in Fig. 2] at a single point on the con- 
trol landscape, it is worth mentioning that the correlator 
identifies such regions on the landscape with remarkable 
ease. 

A noteworthy feature on both the control landscapes is 
the existence of "walls" of no control. Interestingly, the 
wall at Ai « 9.2 is particularly robust and, hence, consis- 
tent with the results in Fig. [1] and Fig. 0] which suggest 
extensive participation by delocalized Floquet states. 
Therefore our observations establish C Z (A) to be a highly 
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sensitive measure capable of identifying the lack of con- 
trol due to the involvement of chaotic states. Several 
other regions of no control can be seen in Fig. [6]A_ which 
are less robust in comparison to the one at Ai = 9.2. 
For instance, the walls around Ai » 7.0, and 8.2 seem to 
break up at certain values of A2 but, interestingly, reap- 
pear at higher A2. In every instance where C Z (A) < we 
have confirmed the presence of atleast one chaotic state 
with substantial p zn . This points to a fairly nontrivial 
role played by the chaotic states in the control process. 
The walls and islands of no control result in the highly 
convoluted topology of the control landscape. Hence, us- 
ing the symmetry breaking property of the 2w-field for 
control purposes is not very effective when chaotic Flo- 
quet states are participating in the dynamics. On the 
other hand regions with C 2 (A) > in Fig. [SK indi- 
cate varying degree of control. Interestingly, an earlier 
analysis^, corresponding to case (C) in Fig. [IJ happens 
to be in a region away from the walls of no control! Again, 
the large positive value of the correlator indicates the ab- 
sence of interactions with chaotic states leading to local- 
ization due to symmetry breaking as seen in Fig. |4p. 

Is it possible to break the wall of no control with a 
judicious choice of <p? In order to explore such a possi- 
bility we show the control landscape for different choices 
of <j) in Fig. [BJ3, C, and D. The case <fi = ir/2 is shown in 
Fig. [6)3 and one observes that C Z (X) < over most of 
the landscape which indicates very little control. This is 
consistent with the fact that the combined driving field 
is symmetric about the time axis for <f> = tt/2. The con- 
trol landscape for <p = 7r/3, and 7r/4 are also shown in 
Fig. and B respectively. Note that decreasing cj> from 
the symmetric value of tt/2 gradually leads to regions of 
control. Further confirmation comes from the survival 
probability data shown in Fig. [6] for two points on each 
of the control landscapes with cf) =^ 0. It is intriguing to 
see that the wall of no control seen in Fig.[6K for Ai ~ 9.2 
is present for all the values of 4> shown here. The reasons 
for such robustness is not clear at the moment. Note 
that such regions in Fig. [6] for 0^0 are characterized by 
large negative values of C Z (A). This seems to contradict 
the earlier work o 39 ' 40 which associated a large value for 
\C Z \ with extensive localization in phase space. However, 
as stated earlier, in the present case the interpretations 
are subtle with the sign of the chaotic level- velocity play- 
ing a key role. Our preliminary studies indicate that in 
such regions the assumption (v) r <C {v 2 )V 2 is weakened 
presumably due to an increased fraction of chaotic states. 
Further work is needed to understand this aspect. Never- 
theless, the numerical results presented in Fig.[6]indicate 
that regions of no control are signalled by C z < whereas 
C z > implies control. 



V. CONCLUSIONS 

In summary, the strong field control landscapes in 
Fig. [5] show that nearly perfect control is obtained as 



long as there is no involvement of chaotic states in the 
dynamics. However, involvement of the chaotic states 
leads to lack of control and such regions on the control 
landscape are interspersed among controllable regions in 
a complicated fashion giving rise to a highly convoluted 
landscape topology. Some of these regions are partic- 
ularly robust under field amplitude and relative phase 
variations (cf. Fig. [BJ) . Given that convoluted landscapes 
are expected in OCT studie o 30 i 32 on classically chaotic 
systems, it is reasonable to expect that chaotic states 
can be associated with certain suboptimal features on 
quantum control landscapes. For systems with effective 
one expects 3 ^ several such plateau regions and the 
effectiveness of the control would depend sensitively on 
the initial state and field parameters. 



We emphasize that the association of uncontrollable 
regions on the landscape with influence of the chaotic 
Floquet states is made possible solely due to the highly 
sensitive nature of the intensity- velocity correlator C z (A) 
(cf . Eq. [8]) . The sensitivity of C z (A) to the underlying 
phase space structures is evident from Fig.[5]and makes it 
an ideal candidate to understand the origins of subopti- 
mal structures on the control landscape from a classical- 
quantum correspondence perspective. Moreover, C 2 (A) 
can be obtained with a few Floquet calculations or via the 
parametric equations of motion^ 41 - and hence applicable 
to multidimensional driven systems. A further attractive 
feature, from the point of multidimensional systems, is 
that the correlator is capable of establishing the influ- 
ence of the chaotic states on the control landscape with- 
out the need for visualizing/classifying the states in the 
underlying phase space. 



Although the present work deals explicitly with a tun- 
neling system, many of our observations are valid in 
other contexts as well. In general, controlling the de- 
cay (dissociation, ionization) of a state is expected to be 
complicated if one or several chaotic states participate 
in the dynamics. The complications will be reflected 
in the form of the optimal fields. This does not con- 
tradict the recent work-22 by Takami and Fujisaki since 
they are concerned with controlling the transition be- 
tween two random states in a completely chaotic system. 
It is the mixed regular-chaotic phase space scenario, a 
generic one for molecular systems, which is expected to 
exhibit some of the features seen in the current study. 
Finally, it would be interesting to see if the influence of 
chaotic states are indeed reflected in the topology of the 
quantum control level sets in an OCT calculation. In 
this regard the present work suggests that the intensity- 
velocity correlator could prove to be a useful tool. We 
are currently exploring the usefulness of the correlator in 
other systems and as a tool in formulating local control 
strategies. 
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